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Abstract 

We develop a theory for a novel state of 4 He films, that possesses off diagonal 
order (as in the superffuid state), as well as hexatic or bond orientational 
order. Within our description, both the hexatic and superfluid transitions are 
still of the Kosterlitz-Thouless type, but the superfluid stiffness is sensitive to 
the hexatic transition. We briefly discuss the possible relevance of this work 
to recent experiments on submonolayer helium films. 
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In the early 1970's a substantial amount of theoretical work was done on the topic of 
supersolid helium. [|I|-|T0|| Such a system would possess both the quantum coherence of the 



superfluid state, characterized by off-diagonal long-range order in the density matrix (or 
"ODLRO" |11|), and crystalline order, characterized by long-range order in the density at 



some finite wavevector, or "diagonal" long range order in the density matrix. Progress in 
the field was hampered by a lack of convincing experimental evidence that such a state 
existed in nature. In this Letter we investigate the properties of two-dimensional helium 
films which have off-diagonal order as well as hexatic order (characterized by algebraic long- 
range order in the orientation of directions to nearest-neighbor atoms). We call such a 
system a "superhexatic," to distinguish it from the superfluid state which possesses no such 
orientational order. 

Superhexatic systems are interesting for a number of reasons. First, the hexatic con- 
straints required on the many-body wavefunction are prima facie much weaker than the 
crystalline ones of a supersolid, @ which should make it easier for the quantum and spatial 
order to coexist. Second, superhexatics can be mapped to a variety of other systems of 
physical interest such as the two-dimensional XY magnet on a random lattice, vortex lines 
in superconductors, or a set of two-dimensional Coulomb charges in quantum gravity. 
|12|1 In such models there are two different species of charge (for example, point masses and 



Coulomb charges), each with its own dynamics and each influencing the other. Third, the 
helium system not only provides an additional theoretical tack, but there is some evidence 
that superhexatics may have been seen in experiment. Recent third sound measurements 
of submonolayer films of helium on hydrogen and deuterium substrates demonstrate two 
independent Kosterlitz-Thouless (KT) transitions: the standard superfluid transition near 
IK and a second one near 0.5K. fT3| , [14|j As we shall discuss below, the lower transition may 
be the freezing of the superfluid to form a superhexatic as the temperature is lowered. 

We shall show below that couplings between off-diagonal and hexatic order can exist, 
and that they need not destroy the superfluid or hexatic transitions. These couplings do 
allow for evidence of the hexatic transition to be seen in the superfluid stiffness. We do not 



address the microscopic interactions that give rise to these couplings in this work. Finally, 
we make a brief comparison of this work to experiment. 

A qualitative [|T5|-|i8| and quantitative [[[JJ theory of the superfluid phase transition can 
be built up from the theory of "ring exchanges" . In this approach the helium atoms possess 
only their kinetic energy and a short-range repulsive interaction. The partition function 
is calculated using standard path-integral techniques, in which the system is evolved for 
an imaginary "time" Ti(3 = h/kT, where T is the temperature. The indistinguishability of 
the bosons allows for contributions to the partition function in which helium atoms have 
exchanged positions. The superfluid transition occurs when it becomes entropically favorable 
for large "rings" of bosons to permute their positions, establishing phase coherence across 



the system. [17 



Hexatic order can affect this process in at least two ways. First, the additional stiffness 
of the system might reduce the amplitude for helium atoms to exchange position, p|j20 



thereby drastically lowering the superfluid density. Second, there is a subtle topological 
effect. The hexatic-fluid transition can be viewed as a disclination-unbinding transition, 
[2T|-p3] where a perfect disclination is a topological defect in the bond-orientational order. 



These point disclinations can interact with the point charges in the superfluid order. The 
vortex-disclination interaction appears most naturally in the context of an XY model in 
a fluctuating geometry. Consider, by way of example, a lattice of helium containing a 
single disclination. The theory of ring-exchanges of the atoms can be mapped to a Landau- 
Ginzburg theory for the superfluid. The energy depends only upon gradients of the phase, 
and can be related to two-dimensional electrostatics |22| , |24|| wherein vortices play the role 
of Coulomb charges. The effect of the disclination can then be viewed as distorting the 
plane containing the charges into either a cone for a negative disclination, or a saddle 



for a positive disclination. [^] Using conformal mapping one can solve the electrostatics 
problem for a single charge on a cone (saddle) and find that it is repelled (attracted) to 
the disclination independent of the sign of the charge, and that the energy of interaction 
depends logarithmically in their separation. A similar calculation shows that the energy of 



a vortex and a dislocation varies inversely with the distance between them. 

It is also possible to demonstrate couplings between the hexatic and superfluid order using 
a microscopic many-body approach, [^6j but here we will consider only a phenomenological 
Landau- Ginzburg model. The free energy of a superhexatic may be written as a sum of the 
elastic energy, the superffuid energy, and an interaction energy. The elastic term is given by 
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where \i and v are elastic constants, xif) is the "gauge field" of the stress tensor, a^, so 
that (Tijif) = eikejidkdix(r) and 6(f) and 0(f) are the dislocation (Burger's vector) and 
disclination densities with core energies E h and Eq, respectively. f27j These defect densities 
can be combined into a single scalar, the incompatibility, 77(f) = 0(f) +z- V x 6(f) which acts 
as a source term for the field x(f) , where z is a unit vector normal to the plane of the film. We 
assume that we are well above the melting transition of the solid, and treat the dislocation 



density as a continuous field, [^TJ ignoring the discrete nature of the dislocations. We next 
write the superfluid energy in terms of a ficticious electrostatic potential 0(f), related to the 
phase of the superfluid order parameter, 6(r), by V0(f) = —z x V0(f). Using this potential, 
the superfluid energy may be written in a fashion similar to the elastic energy: 
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where v(r) is the density of point vortex "charges" for the field ip(r). 

The leading order interaction between the two fields may be written as: 



Etot = J dr il'jkemdjdkX(r)d t ip(r)d m tp(f) + A j H m 6 j (f)6 fc (f)9^(f)9 m ^(f) + E Qv Q{f)v (f). (3) 

The coupling tensors 7^ and \ju are required to be symmetric in order to preserve rota- 
tional invariance. They may be split into the traceful and traceless contributions. For the 
purposes of this calculation we consider only the trace (diagonal) contributions, 
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This is the lowest order interaction one can write for the two fields. Time reversal invariance 
places a powerful constraint on the theory: any interaction requires an even number of powers 
of Vyj; there is no similar requirement on V%(f). The first term in eq.(|]) represents the 
coupling of the stress tensor to the superflow. In the absence of dislocations it generates the 
logarithmic interactions between vortices and disclinations mentioned earlier. The second 
term represents the possible variation of the local superffuid density (stiffness) with the 
local dislocation density. The final term represents the interaction between the cores of 
disclinations and vortices, and can be shown to be generated from the first term in eq.|j. 
The partition function of the superhexatic system may be written as: 

Z= [V(p[r\V X [r\Vb[r\e- f3 ^ +E ^ +E ^ ) (5) 

We simplify the partition function in two steps. First, we integrate out the gaussian dis- 
location vector field, b(f), treating Ao as a small parameter in the resulting action. This 
introduces a |Vx(j)| 2 term in the action so that the interaction of bare disclinations is 
reduced from r 2 logr in the absence of dislocations, to a logarithmic one due to partial 



screening by the dislocations. |2T] Second, we will work in the small fugacity limit, so that 
we may limit the contributions of the sums over vortex and disclination charge to and ±1, 



turning the problem into that of two coupled sine-Gordon systems, pj] After some algebra 
we obtain: 
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The constants gi can be simply related to the original core energies: g\ and g 2 are the 
fugacities of the vortices and disclinations, g 3 represents the added energy required to place 
a vortex on a disclination, and g^ reflects whether or not the energy cost depends upon the 
sign of the disclination. 

Let Te and T v be the hexatic and superffuid transition temperatures, respectively, for the 
uncoupled models (T u > T@). Elementary power counting indicates that all of the couplings 



between the two order parameters allowed by time-reversal symmetry are irrelevant (in the 
renormalizat ion-group sense) at the gaussian fixed point, Tq. A more rigorous calculation 



using momentum shell cutoff renormalization |28| upholds this conclusion near the line of 
gaussian fixed points, even in the regime Tq < T < T u , where the disclination fugacity g 2 is 
perturbatively relevant. |3B| This result is quite important: it shows that weak coupling of 
the two models will not destroy the KT nature of either the vortex-unbinding or disclination- 
unbinding transition. It also allows us to conclude that there will not be a discontinuous 
jump in the superfluid stiffness when the hexatic transition occurs. 

However, this does not prove that the superfluid stiffness is wholly insensitive to the 
presence of the hexatic order. For example, when we pass through the hexatic transition we 
expect a change in ^| Vx| 2 ) proportional to the integral of the bump in the specific heat. Such 
a bump is a non-universal feature of the transition in that its width and magnitude depend 
upon detailed features of the system. p9[ If we have a coupling of the form |Vx| 2 |V<^| 2 , 



then we would expect a change in the superfluid stiffness proportional to this integrated 
bump. For a suitable choice of parameters this can be fairly sharp. However, the lowest- 
order perturbative renormalization group analysis discussed above is only valid near the 
transition, and cannot demonstrate such a feature. 

As a proof-of-principle, we have simulated the superhexatic transition numerically using 
a Monte Carlo analysis. We can define a hexatic order parameter ijj§(r) = ^neighbors e%Wl = 
| ^6 (01 e 106 ^, so that 6* 6 (f) plays a role similar to that of the superfluid phase 9(f), and like 
the phase, its gradient is at right angles to the gradient of its corresponding gauge field, 
Vx(0> introduced above. We simulate the problem as two coupled XY models on a discrete 
lattice: 

E = ~ E fc i cos - + h cos (0 6 ( n ) - 6 (tv)) 

(id) 

+s [cos(0(ri) - 9( rj )) - 1] [cos^efc) - 6 fa)) - 1] (7) 

For a suitable choice of parameters we can obtain a rapid change in the superfluid stiffness 
as a function of temperature near Tq, as shown in fig.(l). 

6 



We conclude with a brief comparison of this work with recent experiments by Cheng and 
Mochel |l"3| , [T4|j as well as other theoretical approaches to the same experiments. ||30|J3T[1 Mea- 



surements of third sound velocities in submonolayer 4 He films on hydrogen and deuterium 
substrates suggest the existence of two KT transitions. The upper transition occurs near 
IK, and is consistent with the standard superfluid transition on conventional substrates. 
The lower one occurs at roughly 0.5K, and is signalled by a change in the third sound veloc- 
ity as the temperature is lowered: for deuterium substrates the velocity drops, whereas for 
hydrogen the velocity increases. The temperature of this transition scales linearly with the 
4 He density. For films thicker than a single monolayer deposited on hydrogen substrates, an 
additional third sound mode appears over the temperature range between the two transi- 
tions. 

If the lower transition is indeed KT, then it seems likely to be a melting transition of 
either a solid or a hexatic. A superhexatic would have no shear modulus, so that superflow 
would not be reduced by pinning, which is consistent with experiment. Depending upon 
the details of the coupling between the hexatic and the superfluid order parameters, it is 
possible to obtain an increase or decrease in the superfluid stiffness. This change may be 
sharp, but is not discontinuous, which is also consistent with the experiments. 

The superhexatic model cannot explain the additional third sound mode seen in some 
cases. However, this mode is only seen for coverages greater than one monolayer, and may 
result from different dynamics in the two layers. A more crucial difficulty is dealing with the 
effect of the substrate potential on the hexatic transition. In classical systems a hexagonal 
substrate will lock in the hexatic order, so that no transition occurs. The large zero 



point motion of both the helium and the hydrogen may mitigate this effect, but that must 
be demonstrated by further study. 



Instead of considering spatial ordering of the helium atoms, Zhang []30|] and Kapitulnik 
HJ have independently postulated an ordering of the thermally excited vortices. In such 
a picture the intermediate state would be a vortex-antivortex lattice (VVL) with sponta- 
neously broken time-reversal symmetry. In this picture the upper transition is the melting of 



the lattice to produce a normal fluid, and the lower transition is the sublimation of the lattice 
as its constituents disappear. Zhang has shown how the additional third sound branch can 
be nicely explained in terms of the optical modes of the VVL. However, in such a picture 
neither transition should produce KT behavior for the superfluid. The upper melting transi- 
tion would produce a universal jump in the shear modulus of the VVL, not in the superfluid 
density, as is observed experimentally, and as predicted in our model. In addition, there is 
no reason why the additional branch should only be observed on hydrogen subtrates when 
the coverage exceeds one monolayer. 

The two pictures differ in that the superhexatic phase should occur below the new, lower 
transition, while the VVL is predicted to occur at and above it. By measuring the vortex 
diffusivity and the onset of nonlinear dissipation as a function of temperature, one should 
be able to establish which regime exhibits anomalous ordering. 

The above analysis of the interaction of spatial and superfluid order can be extended to 
physical systems with a similar mathematical description. In particular, it can be applied 
to an array of Josephson junctions that has been deliberately constructed to include a 
disclination or dislocation. The geometric analysis discussed above would predict a force on 
vortices that attracts or repels them from the defect center. Such an interaction might be 
visible in experiments on the ballistic motion of vortices. 

In conclusion, we have introduced the notion of a novel phase for 4 He films called a 
"superhexatic", which displays both off-diagonal and bond-orientational order. We have 
shown that couplings between the two exist, and that the couplings preserve the KT nature 
of the two transitions. We have shown that this phase is roughly consistent with experiments, 
and suggested ways to test if this approach is correct. Study of the system dynamics is 
underway. [[26|] 
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FIGURES 

FIG. 1. The jump in the superfluid stiffness during the hexatic transition, as determined 
by a Monte Carlo simulation of the action in eq.(^) for a 100x100 lattice with periodic boundary 
conditions. The size and width of the jump are proportional to that of the integral of the specific 
heat of the hexatic over the transition. The parameters K\ = 1 and K2 = 2, were chosen so that 
the uncoupled hexatic and superfluid transitions occur at 1 and 2, in these dimensionless units. 
The stiffness is calculated for coupling S = (circles) and S = 0.5 (squares). 
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